Practical design and simulation of silicon-based quantum dot qubits 
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Spins based in silicon provide one of the most promising architectures for quantum computing. 
A scalable design for silicon- germanium quantum dot qubits is presented. The design incorporates 
vertical and lateral tunneling. Simulations of a four-qubit array suggest that the design will enable 
single electron occupation of each dot of a many-dot array. Performing two-qubit operations has a 
negligible effect on other qubits in the array. Simulation results are used to translate error correction 
requirements into specifications for gate- volt age control electronics. This translation is a necessary 
link between error correction theory and device physics. 
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Quantum computation would enable huge speedups 
of certain very hard problems, such as factorization. 1 
However, quantum computing is essentially an analog 
method. As such, the problem of errors creates a se- 
rious challenge. Advances in error correction algorithms 
have produced well-justified optimism that this challenge 
will be overcome. 2 However, existing error correction al- 
gorithms require low error rates. Thus, hardware design 
will be critical to the creation of a working quantum com- 
puter. 

The purpose of this paper is to address hardware de- 
sign challenges in a specific materials system: silicon 
quantum dots. There are two reasons to analyze this 
system in detail. First, spins in silicon have long co- 
herence times^ Second, classical silicon electronics has 
demonstrated fast operation and a proven record of scal- 
able integration. Indeed, several spin-based qubit designs 
have emerged that are compatible with silicon 
The full benefit of existing silicon technology may be 
used to greatest advantage in spin qubits based in quan- 
tum dots. Previous calculations of the exchange coupling 
in coupled quantum dots with idealized potentials have 
demonstrated the promise of such structures for quantum 
computation However, there are important questions 
that can only be addressed in the context of an explicit 
physical design and realistic simulations. 

In this paper, we present an explicit design for quan- 
tum dot qubits in silicon-germanium heterostructures. 
To determine whether it will be possible to build and 
operate such a device, we perform realistic simulations of 
four coupled qubits. The simulations are self-consistent: 
they include the full three-dimensional electrostatics, and 
the Hamiltonian is solved via exact diagonalization^i 
These simulations allow us to answer several questions. 
First, we find that it is possible to couple neighbor- 
ing quantum dot qubits without any significant pertur- 
bation of secondary qubits. Second, the coupling can 
be strong, enabling GHz operation rates. Most impor- 
tantly, these simulations allow us to translate gate volt- 



age uncertainties — which are inevitable — into error rates 
in quantum gates. This translation is the necessary link 
between device physics and quantum error correction the- 
ory. 

In this paper, we do not propose a new scheme for 
quantum computation. Rather, we perform simulations 
of a new design suitable for implementing the scheme of 
Loss and DiVincenzo^ The quantum computer we have 
in mind is the following: the physical qubits are indi- 
vidual electron spins in quantum dots. Two-qubit op- 
erations are performed on these physical qubits by con- 
trolling the exchange coupling J as a function of time. 
Logical qubits can be coded into a subspace of the phys- 
ical qubits, so that the exchange coupling alone enables 
universal quantum computation™^ Initialization of the 
coded qubits is performed according to the scheme of Di- 
Vincenzo et alm& Readout is performed via spin-charge 
transduction, as in the tunneling scheme of Kane. 5 

The quantum computer just described is well defined, 
but it is abstract. Our specific design is shown in 
Fig. 1(a). It incorporates aspects of two existing types 
quantum dots; lateral tunneling dots and vertical tunnel- 
ing dots. 14 The quantum dot of Fig. 1(a) is defined by 
a quantum well that confines electrons vertically, and by 
split top gates that confine electrons laterally, by elec- 
trostatic repulsion. These features are typical of lateral 
quantum dots. The device of Fig. 1(a) differs from a 
typical lateral quantum dot because it contains a tunnel- 
coupled back gate, usually found only in vertical quan- 
tum dots. As the simulations below show, the back gate 
allows tuning of the electron number in each quantum 
dot, even when those quantum dots are part of a large ar- 
ray. The back gate also screens the Coulomb interaction. 
All semiconductor layers in this design are composed of 
strain-relaxed Sii-^Ge^ except the quantum well, which 
is formed of strained silicon. Relaxation in SiGe can be 
achieved by step-graded compositional growth on a sub- 
strate silicon wafer. 15 In the simulations presented here, 
we use the composition x = 0.077, consistent with a 
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FIG. 1: The two quantum dot devices simulated in this pa- 
per, (a) A double-dot structure, as studied in Figs. 2 and 3. 
From bottom to top, the heterostructure cross section is com- 
posed of a thick, n-doped, strain-relaxed Sii-xGex back-gate, 
a 10 nm undoped Sii-a^Ge^ tunnel barrier, a 6 nm undoped 
Si quantum well, a 20 nm undoped Sii-ccGe x barrier, and 
lithographically-patterned metallic top gates. All fabrication 
steps are based on standard technology. Not pictured is a 
thin Si capping layer, (b) A four-dot structure, as studied 
in Fig. 4. [Top- view only; heterostructure identical to (a).] 
Periodic boundary conditions are assumed. The dots reside 
in the quantum well layer, at positions marked by x 's. 



quantum well band offset of AE C ~ 84 meV with respect 
to the barriers. Since the transverse effective mass and 
dielectric constant change little with x when x is small, 
we use constant values m t = 0.1 9m e and e = 11.9£o 
throughout the heterostructure. 

The zeroth order requirement for an individual elec- 
tron qubit is that it should contain an individual elec- 
tron. Fig. 2 shows the stability energy (the energy cost 
to change the electron number) calculated for the pair 
of coupled quantum dots shown in Fig. 1(a). The re- 
sults are plotted as a function of the gate voltages V\ n 
and lout, as explained in Fig. 3 (inset). The stability 
energy is greater than 1 meV over a wide range of gate 
voltages. We have also performed stability calculations 
for the four-qubit device of Fig. 1(b), and we find stabil- 
ity greater than 1 meV throughout the operating range 
discussed below. 

Qubits are useful only if operations can be performed 
on them. For the structure we describe here, the opera- 
tions are performed by controlling the exchange coupling 
J(t) between neighbor qubit pairs, with the interaction 
Hamiltonian H s (t) = J(t)Si • S2. The exchange coupling 
is large only when the electron wavefunctions overlap. It 
can be made exponentially small by forcing the electrons 
to separate. These manipulations are performed via the 
top gate voltages, and these gate voltages translate di- 
rectly into the time evolution of the qubits. The map- 
ping J(Vi, V2, • . •) between the exchange coupling and the 
top-gate voltages is an operational characterization of the 
quantum computer. As we show below, knowledge of this 



mapping allows us to determine error rates for our design. 
Here, we calculate the exchange coupling as the energy 
difference between the ground and first excited states Aifi 
J = E tr i p — £? S ing, where "singlet" and "triplet" refer 
to the spin symmetry of the two-electron wavefunction. 
Significant numerical accuracy is required in the calcu- 
lations, because of the large difference in energy scales: 

Wrip < 5 X 10- 4 . 

Figure 3 shows a map of the exchange coupling J as 
a function of gate voltages V in and lout, computed nu- 
merically for the double-dot device of Fig 1(a). The 
back-gate is grounded. The envelope function approxi- 
mation used here is reasonable for quantum dots of size 
~ 50 nm£ The overall trends in Fig. 3 are consistent 
with previous studies, which use more idealized confine- 
ment potentials AiS However, because the magnetic field 
is zero in this work, the exchange coupling does not cross 
zero, in contrast with results for high magnetic fields. 
Nonetheless, J can be made arbitrarily small by raising 
the electrostatic barrier between the qubits. Raising such 
a barrier creates an asymptotic approach to zero that is 
extremely robust. 

The data in Figures 2 and 3 are for two qubits. It is 
conceivable that adding additional qubits would cause at 
least one qubit to become unstable or suffer an undesired 
coupling with a neighboring qubit during the manipula- 
tions described by Fig. 3. Fortunately, this is not the 
case. Fig. 4 shows the electron density in a four-qubit 
device for two extreme cases, in which the exchange cou- 
pling is either (a) very small or (b) very large. Between 
Fig. 4(a) and (b), the inner pair of electrons each move 
by 21.5 nm, whereas the outer electrons move by only 0.5 
nm. This motion corresponds to a change in J for the in- 
ner pair from approximately 10 -19 eV to 0.4 /xeV. We can 
estimate the J coupling between the fourth and "fifth" 
electrons in Fig. 4 (using periodic boundary conditions). 
We obtain approximately 10 -19 eV for Fig. 4(a), and this 
number decreases by only a factor of 0.8 in Fig. 4(b). 




FIG. 2: The stability energy (energy to change the electron 
filling number in a two-qubit device) vs. gate voltages, com- 
puted for the two-qubit device, Fig. 1(a). The two-electron 
stability range is shown in the center, while the one-electron 
stability range is shown on the right. Three-electron stabil- 
ity does not occur in the voltage range shown here. Optimal 
two-electron stability is obtained along curve AB. 
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FIG. 3: Exchange coupling J vs. gate voltages, computed for 
the two-qubit device, Fig. 1(a). Top gate voltages Vout and 
Vin are described in the inset, while the back gate is held to 
ground. Curve AB marks the line of maximum stability for 
gate operation. (See Fig. 2.) 



Thus, any pair of qubits can be manipulated indepen- 
dently of any other pair. This independence is due, in 
part, to screening effects arising from the various gates. 

The results of Fig. 3 allow us to consider errors in 
quantum gates. It is important to remember that errors 
arise in quantum computing not just from decoherence 
but also from the inevitable misapplication of quantum 
gates. Such misapplications will arise, for example, from 
uncertainties in the applied gate voltages. Fault-tolerant 
techniques have been developed for correcting errors, but 
these are only effective below an error threshold of one 
accumulated error in 10 4 operations. 17 Thus, it is critical 
to know the error rate expected during the application 
of quantum gates. 

Here we calculate the error rate in J as a function 
of the uncertainty or noise in the voltage pulses used 
to manipulate the quantum dots. Accurate gate control 
involves two steps: (i) initial characterization of the ex- 
change coupling between pairs of qubits, and (ii) precise 
implementation of the gate operations. For this discus- 
sion, we assume perfect characterization, and we focus on 
step (ii). As a prototype for gate operations we consider 
a/SWAP, as implemented with a voltage pulse V s (t). In 
principle, the particular shape of V s (t) is arbitrary, al- 
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FIG. 4: Simulations of the four-qubit device, Fig. 1(b) (gates 
outlined in white), (a) Coupling "off," with J ~ 10 
(top-gate voltages all set to -0.15 V). (b) Coupling "on 
J — 0.4 fieV (center gate voltage set to -0.075 V). 
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though it must satisfy the following relation: 4 

J J[V s (t)]dt = 7Th/2. 



Here, r s is the switching time, and the function J(V) was 
computed in Fig. 3. To be specific, we consider low and 
flat voltage pulses, such that errors in the pulse width 
are diluted to acceptable levels™ 

What are the error levels that can be tolerated in 
the applied gate voltages? For a flattop pulse of height 
J = 7t^/2t s , fault tolerant computation requires that the 
pulse height uncertainty 8V should satisfy 



SV < 10 J 



dJ 



dV 



(2) 



(1) 



It is important to note that the magnitude of 5V depends 
on classical control electronics, while J and dJ/dV are 
implementation-specific. To evaluate gating errors, it is 
therefore necessary to work with a realistic device design. 
By fitting the exponential dependence of J(V) in Fig. 3, 
Eq. (0) yields the requirement SV/V < 5-8 x 10 -6 for the 
double qubit device of Fig. 1(a). We can compare this 
figure to published specifications for state-of-the-art con- 
trol electronics. For sub-kHz pulses (approaching DC), 
extremely high voltage accuracy can be achieved, and 
the requirement can be met. For sub-MHz pulse gener- 
ators, the desired accuracy levels fall nearly within the 
specifications of off-the-shelf electronics™ For GHz op- 
eration, over three orders of magnitude improvement in 
pulse height uncertainties will be required to meet the re- 
quirements of fault tolerant computation™ We point out 
that decoherence constraints may indeed require spin- 
based silicon qubits to operate in the GHz regime. 

In conclusion, we have described and simulated a re- 
alizable design for a SiGe quantum dot quantum com- 
puter. A prominent feature of this device is the back- 
gate, which enables tuning of the number of electrons in 
each quantum dot. We have directly addressed the issue 
of scalability through simulations of a four-qubit device. 
Qubit interactions are found to be very robust, partic- 
ularly as a consequence of Coulomb screening provided 
by the back-gate. Our calculations show that a key chal- 
lenge for solid-state spin-based quantum computation is 
to develop devices in which the exchange coupling is rela- 
tively insensitive to gate voltage uncertainty. At a simple 
level, the quantum dot structures should be optimized to 
increase J/\dJ/dV\, which sets the scale for gate volt- 
age accuracy requirements. The ultimate goal should be 
to "digitize" the gating function J(V), such that dJ/dV 
goes to zero at appropriate working points. 21 
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